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1. Introduction 

Recently the problem of the behavior of the fields in a cavity with a (periodically) moving boundary has 
received significant attention. From a mathematical point of view, it constitutes an instructive example of 
parametrically driven system that exhibits interesting resonant effects. From point of view of physics, besides 
its fundamental importance as a modification of Casimir effect, it can be used as a model for the mechanism 
of formation of wave packet in lasers, processes in atomic physics, and even interstellar flight [ITj! We would 
like to draw reader's attention to the reviews of Casimir effect by Bordag et al [3], the books by Milton |5(Jj . 
and Mostepanenko and Trunov j^, and to the recent review by Dodonov [201 devoted specifically to the 
dynamical (or nonstationary) Casimir effect (i.e., the quantum effect in a pulsating cavity) which contains 
more than 300 references. 

In this paper we will apply a method we have developed in 45 (and generalized in |54| 1 to study 
the behavior of the classical electromagnetic field in a one-dimensional cavity with a moving wall by using 
methods of dynamical systems. Here we will employ this methodology to analyze the quantum problem. Our 
approach is based on studying the collective behavior of the characteristics of the wave equation by applying 
theory of circle maps. The general theorems allow us to predict the behavior of the system without solving 
partial differential equations. We show that the mechanism of the resonant amplification of the quantum 



field is the same as in the classical case - it is due to Doppler effect at reflection from the moving mirror. 
In the quantum case, however, the motion of the mirror "creates" new field which is then amplified by the 
Doppler effect. Using some mathematical identities, we give simple physical interpretation of the different 
contributions to the energy density. 

In the rest of the Introduction, we review very briefly some of the literature related to the mathematical 
and physical aspects of our approach, referring the reader to the review |2Uj . 

The mathematical theory of the solutions of the wave equation in presence of (periodically) moving 
boundaries in one or more spatial dimensions has been developed by Cooper 0011], Cooper and Koch 
[3], Yamaguchi and collaborators [SHI EOl E21 Ej j Dittrich et al JTj. In the physics literature, the classical 
version of the problem was studied by Dittrich et al ^Hl) Cole and Schieve 0, Meplan and Gignoux who 
used the geometric method of solving the wave equation (i.e., the method of characteristics), and recognized 
that under in some cases the fleld develops wave packets that become narrower with time. 

The foundations of the quantum theory of the problem of a 1-dimensional resonator with a moving 
wall were laid by Moore [HI], and developed by Fulling and Davies EHllIll- Dodonov et al 12312311121 
considered the case of resonant motion of the mirror within Moore's formalism; they and Jaekel and Reynaud 
P^. Meplan and Gignoux and others predicted that the force between the mirrors can be enhanced 
significantly in the resonant case. 

In j42| . Law proposed an exact analytic solution for a particular choice of the motion of the mirror, for 
which the field in the cavity develops two wave packets which become narrower in time and whose energy 
grows. Law's solution was generalized by Ying Wu et al 58 , who constructed motions of the mirror for 
which the field develops several wave packets. 

Dodonov ^ni noticed that small enough "detuning" from the exact resonance conditions does not change 
the qualitative features of the behavior of the field in the cavity. Within our approach, one can find explicitly 
for what detuning the behavior of the field will change dramatically. The packet formation in a resonantly 
pulsating resonator was studied in detail by Dodonov and Andreata |^ ^ . 

Ideas from dynamical systems (i.e., study of the iterates of a certain map) have been used by Cooper 
13 m (maps of the interval) , Meplan and Gignoux (area-preserving maps of the cylinder) , Dittrich et al 
[Hini, Yamaguchi (5112011221 (circle maps). 

2. Method of characteristics and dynamical systems 

In this section, we explain the physical setup, pose the mathematical problem, and explain how to analyze 
it. We recommend that the reader consult our paper [32] for details. 

We will refer often to the electromagnetic (EM) field in the cavity (meaning the classical, not quantum, 
field); as it turns out, some aspects of the problem are similar in the classical and in the quantum problem 
case. 

2.1. Description of the physical system 

Consider the EM field in the empty space (no medium, electric charges, or currents) between two parallel 
perfectly refiecting mirrors, one stationary at a; = 0, the other moving according to a; = a{t). The function 
a must satisfy the physically natural conditions a{t) > for each t (the resonator never collapses to zero 
length), \a'{t)\ < 1 (the speed of the moving mirror never exceeds the speed of light). To avoid technicalities. 
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we assume that a{t) is a smooth (C°°) function. Wc wih focus on the case in which the motion of the mirror 
is 1-periodic (i.e., periodic of period 1): 

a{t + 1) = a{t) for each t . (1) 

Examples of such functions are 

ait) ^ - + ^ sm2TTt , ■7^>^, |/3|<1; 

«W = f + |^sin(27rt + 7(sin4^t)2) , | > M , + aj^l) < 1 . (2) 

In our numerical simulations we will use the function a{t) from (j^J. 

2.2. Boundary-value problem 

Since there are no charges and currents in the cavity, we impose Coulomb gauge = 0, V • A = on the 
EM 4-potential = (Aq, A), and obtain that A satisfies the wave equation. We consider linearly polarized 
plane waves propagating in x direction, so we can assume without loss of generality that the vector potential 
has the form 

A{t, X) = A{t, X) By . 

The function A{t, x) satisfies the (1 + l)-dimensional wave equation 

Att{t,x)~A,,{t,x)^0 (3) 

in the spatio-temporal domain A := {{t, x) G M^|0 < x < a{t), t > 0}. The boundary conditions (BCs) come 
from the fact that in the coordinate frame instantaneously co-moving with the mirror, the tangential to the 
mirror component of the electric field must vanish at the mirror, which yields the "perfect reflection" BCs 

At(i,0) = 0, At{t,a{t)) + a'it)A,{t,ait))^0 (4) 

for each t > 0. Geometrically, the BCs mean that the derivative of A{t,x) along the world line of the 
mirror (i.e., the line {{t,a{t)) | t £ R} in the space-time diagram) must be 0. Note that the Dirichlet BCs 

A{t, 0) — ci — const , A{t, a{t)) = C2 — const (5) 

are equivalent to Parenthetically, we would like to note that Neumann BCs are not Lorenz covariant, so 
they are not physically natural for the case of EM fields. If, nevertheless, one imposes Neumann BCs, the 
predictions of the theory are dramatically different from those of the "perfect reflection" BCs (see |^ 
Section 5.6] or Dittrich et al ^\ Section 4]). 

In the quantum case (Moore |51)). one obtains the boundary value problem consisting of the equation 

the Dirichlet BCs (|3Jl with ci ~ C2 — 0, and initial conditions. 

2.3. Method of characteristics 

In absence of spatial boundaries (i.e., if A = {(t, x) G E^|t > 0}), the solution of the wave equation J^Jl with 
"perfect reflection" Q or Dirichlet BCs, and initial conditions A{0,x) — ipi{x), At{Q,x) ~ tp2{x), x € M, 
is a superposition of waves propagating to the left and to the right: A{t,x) ~ ^'+(a;(j") + '^^(xq"), where 
Xq = x ±t, and ^'^(s) = ^ ["01 (s) ± ij^iis') dx] (c is an arbitrary constant, the same for ^E'^ and ^'^). 
Geometrically, in the space-time diagram, the waves propagate along the characteristics, x ±t = const. 
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Figure 1. Characteristics of the wave equation. 

In the presence of spatial boundaries (stationary or moving), the characteristics are no more straight 
hnes, but are piecewise hnear, each part of them being a straight hne at 45° with respect to the t axis (see 
Figure In order for the BCs Q to be satisfied, the field changes sign at each reflection, so that 

Ait,x) = (-l)^+vI,+ (4) + (-l)^-vI;-(^-) , 

where N± is the number of reflections of the corresponding broken characteristic by a mirror between the 
initial moment t — and the present time t. For a proof that this algorithm works, see 0^1 Section II. B]. 




2.4- The importance of the reflections; Doppler effect 

Loosely speaking, the "density" of the characteristics in the space-time diagram is proportional to the energy 
density of the EM field. Consider two characteristics corresponding to the ends of a narrow wave packet. 
During free propagation (no reflection), the width of the wave packet and the vector potential A do not 
change. At reflection from the moving mirror, however, not only does A change sign, but also the width of the 
wave packet changes. If the width of the wave packet before the reflection was A, and the reflection occurs 
at time t (assume that the wave packet is so narrow that the reflection happens almost instantaneously), 
simple trigonometry shows that after the reflection the width A' of the wave packet is A' = A/D{t), where 

is the Doppler factor at reflection at time t. The term "Doppler factor" comes from the fact that the initial 
(classical) EM energy of the wave packet, [Atitjx)"^ + A.j;{t, x)'^'^ dx, increases by a factor of D{t) at 

reflection from the moving mirror (see |54[ Section 2.4] for a simple proof). Clearly, if the mirror is moving 
inwards (outwards) at the time of reflection, the energy of the wave packet will increase (decrease). If it 
happens that every time a certain group of nearby characteristics (representing a wave packet) is reflected 
from the moving mirror while the mirror is moving inwards, then they are going to get closer together, which 
will lead to squeezing of the wave packet (see Dodonov et al and Jaekel and Reynaud |^) and to 

exponential growth of the energy of the field, as we will see below. 
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Since the motion of the boundary is 1-periodic the position a{t) and the velocity a'(t) of the mirror, 
as well as the Doppler factor D{t) 0, do not depend on the integer part of t, but only on its fractional part, 

{t}:=t-[t\, (7) 

which we will refer to as the phase of the motion of the mirror. To make the phase of t change continuously 
as t increases, we will think of {t} as belonging to a circle of length 1, i.e., to the interval [0, 1] with its ends 
identified. The long-time behavior of the field depends on the asymptotic behavior of the characteristics, 
which in turn can be analyzed by invoking the mathematical theory of circle maps, as explained below. 

2.5. From characteristics to circle maps 

Since characteristics belong to a very simple class of plane curves - namely, piecewise linear at a 45° angle 
with the t axis, - to reconstruct a particular characteristic, it is enough to know only one moment at which 
this characteristic is reflected from, say, the stationary mirror. 

To study the collective behavior of the characteristics, we introduce the time advance map F : M. ^ R 
such that if certain characteristic is reflected from the stationary mirror at time ti, the next reflection from 
the same mirror occurs at time t2 = F{ti) (see Figure^. To derive an expression for F in terms of the 
function a giving the motion of the mirror, we notice that the time Q{ti) between ti and ^2 at which the 
characteristic is reflected from the moving mirror satisfies — a(0(ii)) = ti, which can be written as 

(Id - a)(e(ii)) = ti, therefore 

e=(Id-a)"^. (8) 

On the other hand, F{ti) = Q[ti) + 0(8(^1)) = (Id + a){Q[ti)), thus 

F^ (Id + a) o (Id-a)"i . (9) 

The conditions a{t) > and \a'{t)\ < 1 guarantee the invertibility of (Id ± a) (hence the existence of 8 
and F) as well as the fact that F is strictly increasing and, therefore, invertible. We leave to the reader to 
check that a can be expressed in terms of F as 

a = i(F-Id)o [i(F + Id)]"' . 

The 1-periodicity |^ of a guarantees that F satisfies the property 

F{t+l) = F{t) + l . (10) 

Since only the phase {t} (TJ is physically important, instead of considering the function F : K ^ M, we 
define the function 

f:S'^S':{t}^ {F{{t})} (11) 

that maps the phase {t} at some reflection from the stationary mirror to the phase {-F({i})} the next 
reflection from the same mirror. Here stands for the "circle", i.e., the interval [0, 1] with its ends identified 
(in mathematical notations, this can be written as — R/Z, where Z stands for the integers). The function 
/ is well-defined due to I|1U|I . 

If the first reflection of a particular characteristic from the stationary mirror occurs at time ii, the times 
of the subsequent reflections are -F(ii), F^(ii), _F^(ti), . . ., where 

F" := F o o • • • o F 
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Figure 2. Graph of the functions F (thin solid hne) and / (thick dashed line) corresponding to the function 
a for a = 0.35, /3 = 0.4, 7 = 0.7. 

is the nth iterate of the function F. Since the asymptotic behavior of the characteristics is completely 
determined by the asymptotic behavior of the phases at reflection, the long-time behavior of the system can 
be studied by analyzing the high iterates of /. The branch of mathematics that studies the behavior of 
highly iterated functions is called theory of dynamical systems. Traditionally, the functions that are going 
to be iterated - like F and / - are called maps. In particular, the map / is an example of a circle map 
(CM), i.e., a map from the circle S-^ to itself. Theory of CMs is a prominent part of theory of dynamical 
systems; it was initiated by Poincare in 1880s, and nowadays is a highly developed field of mathematics with 
many physical applications. 

The relationship between the time advance map : M ^ R and the CM / : S*^ ^ S*""^ is shown pictorially 
in Figure |21 Note that although / looks discontinuous in the figure, it is continuous as a function on the 
circle because of the identifications of and 1 shown in the figure with dotted lines. The map F is called 
a lift of /, while / is sometimes called the projection of F. Clearly, F determines / uniquely; on the other 
hand, each CM / has infinitely many lifts that differ by an additive integer constant (in our case, however, 
the lift F is defined uniquely by ©). 

3. Circle maps and wave packet formation 

In this section we collect some facts about the dynamics, i.e., the behavior of the high iterates /", of CMs. 
For more information the reader can consult the introductory expositions in Hasselblatt and Katok |^ 
Chapter 4] or Devaney ^1 Section 1.14], or the more sophisticated treatments in Katok and Hasselblatt 
Chapters 11 and 12], de Melo and van Strien ^1 Chapter I]. Section III of our paper 0^ contains a selection 
of mathematical facts adapted to the problem of the resonator. In Section IX^ we will give interpretation of 
the mathematical results in terms of the asymptotic behavior of the field in the resonator. 

3.1. Circle maps - basic definitions 

By a circle map (CM), we will always mean a smooth (C°°) invertible map of the circle whose inverse is 
also smooth - this is exactly the class of CMs that correspond to motions of the boundary a satisfying the 
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conditions from Section l2.ll For the map / to be invertible, we have to assume that the cavity is not too 
long, or, more concretely, that F{t) — t < 1 (for which it is enough to assume that a{t) < ^ for aU t). This 
condition only helps to avoid clumsy sentences, but is not a restriction of the generality - our ideas can be 
easily applied mutatis mutandis to the case of a longer cavity. 

The most important characteristic of a CM is its rotation number defined as the "average amount of 
rotation" : 

t(/) ^ t{F) := Km (12) 

71 — >CfD TL 

(with the above restriction on the length of the cavity, t(/) G [0, 1)). It can be proved that r(/) always 
exists and does not depend on the value of t in (|12|l . 

An orbit of a point t E is the set {/"(i)}5^o of future (i.e., for n > 0) iterates of t. If for some 
point t* G there exists an integer q such that fit*) = t*, then we say that t* is a periodic point of 
period q (or a q-periodic point) and call the orbit {/"(i*)}J'i^o '^^ ^^^^^ point a q-periodic orbit. 

The simplest example of a CM is the rigid rotation r„ : ^ (where a £ [0, 1)) defined through its 
lift R^, 

Ra{t):=t + a, teM, 

1 (13) 
r,{t) -.^{t + a} , teS^ . 

Clearly, T{r^) = a. The dynamics of is very simple: 

• if is a rational number, i.e., a = p/q for some integers p and q (we will always assume that p and q do 
not have common factors), then after q iterations any point t G returns to its initial position, having 
traversed the circle p times, i.e., each point t G 5*^ is a g-periodic point: 

Rl/qit) = t + P VteR, 

• if cr is not a rational number, then the orbit {r^ait)}ri'=o of point t ^ fills the circle densely, and 
will never return to the initial point t, thus, in this case there are no periodic orbits. 

3.2. Phase locking, Arnol'd tongues, devil's staircase 

Here we will describe in detail the case of a general CM / with a rational rotation number, t(/) — p/q, in 
which case the map / is said to be phase locked {frequency locked, mode locked). 

If t(/) = p/q, then generically / has an attracting q-periodic orbit {ij'^^lj^i, and a repelling g-periodic 
orbit {tj^^}'j^^. "Attracting" means that the orbit of each point t that is not one of the of the repelling 
periodic points {t'^^jj^i tends asymptotically to the attracting periodic orbit {t'f^j'j^i. The repelling 
periodic orbit "repels" the iterates of /; it is an attracting periodic orbit for the inverse map /^^ (which is 
also a CM) . The attracting and repelling periodic orbits of the CM / give rise to attracting and repelling 
characteristics of the wave equation, and to formation of wave packets (see Section l3.3f) . 

A very important for the physics of the problem question is how "generic" the case of phase locking is. 
In Figure |3Ja) we show in the {a, /3)-plane the regions of values of the parameters a and (3 of the motion of 
the mirror that correspond to phase-locking of several rotation numbers t(/). These "phase- locked" regions 
in the (a, /3)-plane are called Arnol'd tongues in honor of Arnol'd who studied them in his famous paper on 
CMs The Arnol'd tongue corresponding to a. p/q phase locking emanates (i.e., starts as /3 — > 0+) from 
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Figure 3. (a) Arnord tongues with rotation numbers 1/4, 1/3, 2/5, 9/20, 3/5, and 2/3 for the CM / 
corresponding to motion of the boundary with 7 = 0.7. (b) Graph of the rotation number t(/) for 
with P = 0.4, 7 = 0.7. 

a = p/q, and becomes thicker as (3 increases. Tongues with large p and q are very thin (see the 9/20 tongue 
in Figure |3t a)). 

Another illustration of the abundance of phase-locking is the graph of the rotation number t(/) versus 
a (all other parameters fixed), shown in Figure |3Ib). It can be proved that this function is continuous, and 
it is locally constant if t(/) is rational, and strictly increasing if t(/) is irrational, i.e., the graph of such a 
function contains infinitely many densely interspersed horizontal pieces, each corresponding to a particular 
type of phase locking. Such a graph is called a devil's staircase. Several of the horizontal pieces in the figure 
are labeled with p/q showing the type of phase locking. The behavior of the widths of the Arnol'd tongues 
as (3 ^ Q'^ is studied by Jonker and Davie . 

It is worth noting that, in some sense, phase locking is more "generic" than the imlocked case. Namely, 
if t(/) is irrational, then there exists an arbitrarily small smooth perturbation of / such that the perturbed 
map is phase locked. On the other hand, if r(/) = p/q^ then the parameter values are either strictly inside 
or on the boundary of the g/p-Arnol'd tongue. If the parameters are strictly inside the tongue, then a 
small enough (but otherwise arbitrary) smooth perturbation will not change its rotation number, i.e., the 
perturbed map will have rotation number p/q. Practically, however, one should not forget that the width of 
the Arnol'd tongues decreases fast when p and q increase or when (3 decreases. 

If the CM / is in p/q phase locking and the parameters are strictly inside the p/q tongue, / has an 
attracting q-periodic orbit and a repelling one. In this case, if t* belongs to the attracting periodic orbit, 
then f'^{t*) = t* and f'{t*) < 1. In the following, we will say that a phase locking is "generic" if the 
parameters of the CM are strictly inside the Arnol'd tongue, in which case there exist an attracting and a 
repelling orbit. If the parameters of the map are on the boundary of the Arnol'd tongue (i.e., at some end 
on the corresponding horizontal piece of the graph in Figure I^Jb)), then there exists a g-periodic orbit which 
is neither attracting, nor repelling; in this case f'{t*) — 1 where t* is a q-periodic point. 

When P reaches the critical value at which the boundary is moving at the speed of light at some moment 
in each period (i.e., if \a'{t)\ = 1 for some t), then the total length of the phase locking intervals (i.e., of the 
horizontal pieces in Figure |3Ib)) becomes equal to 1, which means physically that the probability of phase 
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locking is 1. This mathematical problem is studied numerically by Jensen et al "6b' and Lanford '41', and 
proved rigorously for general CMs by Graczyk and Swi§,tek . 



3.3. Derivative of the circle map, Doppler factor, formation of wave packets 

Now we translate the mathematical facts about the dynamics of the CM / Hll|) and the time advance map 
F @ into asymptotic properties of the field in the cavity. 

If a particular characteristic is reflected by the stationary mirror at time t, then the times of the 
subsequent reflections from the same mirror are given by F'^{t) {n = 1,2,...); the phases iQ of the motion 
of the mirror at these times are given by the iterates of the corresponding CM, /"({i})- If / phase locked 
with rotation number r(/) = p/q, then, generically, there exists an attracting periodic orbit {t^''''}'^]^ which 
attracts the iterates of any point in (except the repelling periodic points) under the map /. Physically, 

fa) 1 

each attracting periodic point tj € S oi the CM / corresponds to an infinite sequence of times of the 

(a) . . 

form n + tj , where n is any integer, such that the characteristics that are reflected from the stationary 

(a) 

mirror at these times attract the nearby characteristics. We will call times of the form n + tj ' "attracting 
g-periodic times" , and the corresponding characteristics "attracting q-periodic characteristics" . The presence 
of q attracting periodic characteristics means physically that the field in the cavity develops (at most) q wave 
packets, whose widths decrease exponentially, and whose energies increase exponentially with time. 

In the case of generic p/q phase locking, the rate at which the characteristics get closer together is 
related to the first derivative of F (or, equivalently, /), which in turn is related (according to ©, ®) 
to the Doppler factor D{Q{t)) at the time Q{t) of the first reflection from the moving mirror after t: 

^ 1 + a'm)) ^ 1 
1-a'm)) D{e{t))- 

Asymptotically, the wave packets are very narrow, so that they are reflected from the moving mirror 

fa) 

practically instantaneously, at times of the form n + t^ ' . The asymptotic "cumulative" Doppler factor 
Vq over a sequence of q consecutive reflections of the packet from the moving mirror (which takes total time 
p according 
reflections: 



'9 

p according to the fact that F'^ (t'^'^ ) = f j^-* + p) is equal to the product of Doppler factors at each of these q 



Vq:=l[D{Q{tf)) 



llF'it 



(a)N 



{F'^yit'^H' ■ (14) 



This formula holds exactly only in the case of classical EM field, when the motion of the mirror only amplifies 
the field through Doppler effect at reflection. In Section IT!^ we will see that in the quantum case the energy 
is not only amplified, but also created by the motion of the mirror, which introduces corrections to the rate 
of change of the energy. 

In the case when the rotation number t(/) is irrational, the field does not develop wave packets, and 
its energy changes with time, but does not have a tendency towards steady grow or decay. 



4. Quantum effects in a periodically pulsating resonator 
4.1. Moore's functional equation 

Moore ^ was the first to consider the problem of quantizing the electromagnetic field in a one-dimensional 
resonator with a moving wall. We leave out all the complications that he had to overcome in the development 
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of a quantization scheme, and focus on one particular aspect of his treatment (adapting his equations to our 
approach). Let the motion of the mirror correspond to time advance map jnj of rotation number a. Moore 
showed that in this case one has to look for an expansion of the field operator A{t^ x) in mode functions 

=e-2^'*i^(*--) -e-2-'=is(t+.) ^ fc = l,2,3,... , 

where the function S : K ^ M satisfies Moore's functional equation 

+ a{t)) = - a{t)) + a , 

which ensures that A}^(t,x) satisfy the Dirichlet BCs (jsj with ci — 02 — 0. This equation can be rewritten 
in terms of the map as E o F{t) — S](i) + ct, or, equivalently, as 

SoF==i?^oS (15) 

(where H13|) is the rigid rotation by a) , and interpreted as the fact that the value of E changes between 
two consecutive reflections from the stationary mirror by a. This implies, in particular, that if for some 
particular value t we know the values of E(t) in the interval [i, F{t)), then we can reconstruct the function 
E for aU t e M by using ((T^ . 

Moore's functional equation is easy to solve numerically. Let us assume that before t = 0, the two 
mirrors were at rest, and at t = the right mirror started moving: 

{- f < 

' ' ,4 , o^ (16) 

f + 1^(1 -e-* ) sin (27rt + 7(sin 47rt)2) , t>0. 

The motion of the boundary for i > is very similar to the one in (0); we have used the factor 1 — e~* 
to smooth out the transition (i.e., to ensure that E', E", and E'" are continuous, the reason for which will 
become clear in Section IT!^ the factor (1 — e~* ) tends to 1 very quickly as t grows, so it does not affect the 
asymptotic behavior of the system. If the mirrors are stationary, it is natural to take the function T,{t) to 
be linear, so for t G [—a,0), we take E(t) = + const (the constant is immaterial since only the derivatives 
of E have physical meaning), and then use (|15|l to find E(t) for t > 0. 
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Now we will apply our knowledge about the dynamics of the CM / - and, hence, about the time 
advance map - to draw conclusions about the asymptotic behavior of the function S, which will allow 
us to make predictions about the long-time behavior of the energy density of the field. Since the case of a 
rational a = p/q is especially interesting because of the occurrence of resonant phenomena (phase locking), 
we focus on this case in the rest of this subsection. In the case of a generic p/q phase locking, the CM / 
has an attracting and a repelling g-periodic orbits. Let t'^i \ . . . , t'f^ be a sequence of q consecutive attractive 
g-periodic times, and let the ordering be such that 

(this implies that the fractional parts of these times satisfy {t^-f'} • • • {^^q^} ^ {^i^''})- If n. is a positive 
integer, then H15|l iterated nq times reads 

Soi^"9(t) = I](t) + np . 
Using that F"«(tf^) = ^''^ + np, we obtain 

S(tf^ = S(tf^) + np . 

This allows us to define a sequence of functions (for ?i = 0, 1, . . .) 

E„ : + np) ^ E(t«) + np) (17) 

Y.n{t) ■■= S(i + np) ~ np , 

and study the behavior of E(t) for very large t by analyzing the behavior of S„(i) (where t G [t\ ,t\ + np)) 
for n — > oo. The graph of E can be assembled from translates of the graphs of Sn as shown in Figure QJa), 
for parameters corresponding to 1/3 phase locking. In Figure0fb), we show the graphs of several S„ for the 
same parameter values. 

In the case of generic p/q phase locking, the times of reflection, t, F{t), F'^{t), . . ., of a particular 
characteristic from the stationary mirror accumulate at the attracting periodic times (of the form n + t'j^^), 
while the values of T, at two consecutive reflections differ by the constant value a. This difference between 
the behavior of the arguments and the values of the function E (cf. IjlSI) ') explains the occurrence of exactly q 
steep parts of the graphs of S„ at t = t^'^'' (j = 1, . . . ,q), and q almost horizontal parts for large n. Physically, 
the steep parts correspond to the times of reflection of the packets of the field; the "widths" of these packets 
are given approximately by (cf. H14I) ') 

Aj") ^ const. [(^^^)'(4^))r = H^. 
4-2. Energy of the quantum field 

Fulling and Davies '28 computed the energy density of the quantum field in a one-dimensional cavity with 
one stationary and one moving mirror using the "point-splitting" method of DeWitt 15 . They found that 
the regularized energy density in the space between the mirrors (i.e., the energy minus an infinite constant) 
is a superposition of the energies of left- and right-propagating disturbances: 

(roo(t,x)),eg = - ^ Mt + X) + - X)] , (18) 



where 

i2 



27r2 
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and is the Schwarzian derivative of the function S, defined as 



j:"'{z) _ 3 

S'(z) 2 



2 



(19) 



The Schwarzian derivative is a remarkable (highly nonlinear!) differential operator, first introduced in 
complex analysis. If </) is a complex analytic function, then vanishing of is a necessary and sufficient 
condition for to be a Mobius (i.e., fractional linear) transformation, M{z) — ffipj, where ad~ be ^ (see, 
e.g., Nehari [531 Chapter V]). The Schwarzian derivative is invariant with respect to a composition with a 
Mobius transformation, ,5^moF = ^f, which follows from the identity 

yaoHiz) = ^g{H{z)) [H'{z)f + ^h{z) ■ (20) 

The Schwarzian derivative appears in many branches of mathematics - dynamical systems (Singer |55| . de 
Melo and van Strien ^| Chapter 1], Graczyk et al EHI)) Lorentzian geometry (Kostant and Sternberg 
|40|. Duval and Guieu |2S], Duval and Ovsienko Singer theory of differential equations (Hille 

[551 Chapter 10]), integrable systems (Burstall et al 0]), among many others. Even more interestingly, 
Schwarzian derivative is widely used as a tool in theory of CMs (Herman Graczyk and Swi§,tek |30j). 

One can use the the property (|15|l and the composition rule (|20|l to predict the long-time behavior of 
the energy density H18|l. To this end, differentiate both sides of S o — Rj^ o E (which is (|15() iterated j 
times) to obtain 

^'<^'<'>' = WW) ■ 

On the other hand, taking the Schwarzian derivative of the same relationship and using (|20|l with G = E, 
H — , we have 

,y^iF^{t)) = j^-y^[,y^{t)-yp,it)] . 

These expressions yield 

<^iFJ{t)) = 2 Mt) - ^F. W] • (21) 

Using this equation, we can compute the energy of the field at an arbitrary space-time point if the function 
E(i) is known for t G [i,F{i)) where t is an arbitrary value; in particular, if for t < the mirrors are at 
rest, we can take T,{t) as in the discussion after Hl()|) (there we used the smoothing factor 1 — e"'" because 

contains third derivatives of E). In the case of classical EM field, the evolution of the energy density is 
similar to H21|) except for the term S^pj{t), which corresponds to the purely quantum effect of creation of 
field by the motion of the mirror. Of course, in absence of electromagnetic field in the cavity at t = 0, the 
classical energy is zero for all i > 0, while the energy of the quantum field is non-zero for t > even if at 
i = it was zero. 

The composition rule H20|) implies 

fe=0 fe=0 i=0 

which allows us to rewrite equation (|21|l as 

[(F^nt)? f^^^F'{F^-^t)r---F'iFHt)r- ^ > 
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4-.3. Physical mechanism of the energy changes 

In this subsection, we give a transparent physical interpretation of the terms in the right-hand side of H22|l . 
Namely, the first term in the right-hand side of H22|) is the initial energy density amplified in the j reflections 
from the moving mirror between t and (t) , while the term with summation index k in the sum corresponds 
to the field created due to the motion of the mirror at time F'' (t) and subsequently amplified at each of the 
following reflections. To prove this, we will use the following result concerning the energy density emitted 
by a single moving mirror. Fulling and Davies |28| proved that if a perfect mirror is moving in vacuum 
according to a; = ci{t), then the regularized energy density to the right of the mirror (i.e., for x > a{t)) is 
given by 

where F is given by lO, and te is the time of emission. To the left of the mirror (for x < a{t)), the energy 
density is given by 

(roo(t,x)),eg - ^^y^{t + x) = -±-y^{{ld + a){Q) , (23) 

where F ~ (Id — a) o (Id + a)^^ is defined similarly to F 0, but with a replaced by —a. Similarly to ((HJ, 
we define the function O = (Id + a)^^ . 

To rewrite (|23|l in another form, we will need the following formulae: 

(which follows directly from H19|) ). and 

yG{G-\0) = -[G'(G-i(or ^G-i(0 

(a consequence of (1201) )j which in turn implies 

^5((Id + a)(ie)) - -e'((Id + a)(te))VM+a(te) = ~ ■ 

Now we have 

^^((Id + a){Q) = ^(M-a)oe((Id + a)(te)) 

= [e'((Id + a)(te))] Vid-a(ie) + ^e((W + 

= [e'((Id-f a)(te))]' [S^M-a{te) - yid+a{te)] 

^^ a"'{Q[l^a'it,)']+3a'it,)a"{t,r 
[l-a'(ie)]' [l + «'(ie)]' 

Now we will use (|24|l to understand the physical meaning of (|22|l . For simplicity, we consider the 
evolution of the energy density along a characteristic which at time to passes through the point xq while 
moving to the left. At time t+ :— Iq + xq this characteristic is reflected from the stationary mirror, then at 
time &{t+) it is reflected from the moving mirror; the next three reflections occur at times F(t+), 8(F(t_|_)), 
and F^(t_|_), respectively. Let ti be time after F^(t+) but before the next reflection, and xi = ti — F'^{t+) 
be the spatial coordinate of the characteristic at time ti. For the energy density we obtain from H22I) 

(roo(ti, xi)}rcg = D{Q{F(t+))f D{Q[t+)f (roo(to, xo))rcg 

13 



The first term in tlie riglit-liand side of (|25|) is the initial energy density (Too(ioi 2;o)}rcg amplified by 
the factor of D{Q{F{t+)))'^ ~ [F'{t+)]^'^ at the reflection from the moving mirror at time 0(t+), and by the 
factor of D(9(F(t+)))^ — [F' {F{t+))]^^ at the next reflection from the moving mirror at time Q{F{t+)). 

The second term in the right-hand side of (|25|l is the energy density created at the reflection of 
the characteristic from the moving mirror at time and consequently amplified by the factor of 

D(Q{F{t^)))^ at the next reflection from the moving mirror. Indeed, we have 

S^F{t+) = y(M+a)oe{t+) = [e'(i+)]2^id+a(e(t+)) +^e(t+) 
= [e'(i+)]2 [^id+a(e(t+)) -^id-a(e(t+))] 

a'"{0 [l-a'm +3a'{0a"{0^ 



Together with (|24ll . this equality implies 

1 ypiU) _ 1 a"'(0[l-a'(0'] +3a'(0a"(e)' 



«=e(t+) 



24^ F'(t+)2 12^ [1 ^ a'iOr [1 + a'm 

1 



247r 



?=e(t+) 

^^((Id + a)(e(i+))) , (26) 



and a comparison with l|26(l proves the correctness of our interpretation of the term 'F^(t*^)2 ■ 

Similarly, the third term in the right-hand side of l|25|) is the energy density created by the moving 
mirror at time 0{F{t^)). 

This discussion elucidates the difference between the classical and the quantum cases - in the classical 
case the moving mirror amplifies the wave packets by squeezing them, while in the quantum case the moving 
mirror not only amplifies the already existing wave packets, but also creates new field which subsequently is 
amplified at each reflection from the moving mirror. 

4- 4- Resonant amplification in a periodically pulsating cavity 

Now we will consider the particular case of periodic motion of the mirror when the corresponding CM / 
has rational rotation number, t(/) = a — p/q, i.e., is in p/q phase locking. In this case the classical EM 
field in the cavity develops wave packets whose number can be anywhere between 1 and q depending on the 
initial conditions. The wave packets become narrower at each reflection, and their energy increases at each 
reflection. 

In the quantum case, the motion of the mirror itself creates energy which is subsequently concentrated 
in narrow wave packets. If the acceleration of the the mirror is not zero except at isolated times, then the 
mirror is emitting energy all the time, and the number of wave packets developed is exactly q. 

In the case of p/q phase locking, the CM / generically has an attracting g-periodic orbit, which 
corresponds to the times of reflection of the attracting characteristics. If t* is such a time (i.e., if the 
fractional part of t* belongs to the attracting g-periodic orbit of /), then 

F''{t*) = t* +p , < (F«)'(r) < 1 , 

which implies that asymptotically the cumulative Doppler factor Vq ~ [(-F"^)'(i*)]^^ (see ()14|ll is greater 
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than 1. In the case of p/q phase locking, H22|) imphes (for any integer n) 

<^iF-m) - - E -^--^"^"^^^ , . (27) 

If t* is a time of reflection of an attracting characteristic, (|27|) yields 



<f{t* + np) = cO(r ) - Yf, (t*) E 



J=0 



= U(^*) _ <^^,(i*) j . (28) 

The applicability of this expression is not restricted to times like t* - since asymptotically all characteristics 
are very close to the attracting ones, ()28|l gives approximately the asymptotic behavior of the energy density 
of the wave packets. Parenthetically, we remark that using this expression, one can prove that for motion 
of the mirror with parameters corresponding to the ends of the phase-locking intervals (i.e., the ends of the 
horizontal parts of the graph in Figure|3Jb)), the energy density grows not exponentially, but polynomially: 

lim $(t* + np) = <PU*)- nypAt*) . 

The energy of the jth wave packet at time t is given by 

Sj{t) - / (Too(t, a;))rcg dx = --^ / [$(t -x) + + x)] dx , 

where Aj(t) is the support of the jth wave packet at time t. Recalling that at reflection at time t the width 
of the wave packet decreases D{t) times, we obtain that asymptotically the energy of the field in the cavity 
increases exponentially: for large i, 

ftotai(i) ~ const • V'JP . 

The energy density of the wave packets changes after each reflection from the moving mirror, but after 
appropriate rescaling, the its "shape" at times np as n oo tends to a some constant profile which depends 
on the motion of the mirror. In Figure|Sl we show the evolution of the "shape" of a wave packet for mirror's 
motion same as in Figure^a). In the figure we show the rescaled energy density, • (Too(n, a;))i.og, on the 

vertical axis versus the shifted and rescaled spatial coordinate 2?^ • (a: — a;* ) on the horizontal axis, at times 
n for several values of n; here cc* is the spatial coordinate of the attracting characteristic corresponding to 
this packet, at time n. 



5. Concluding remarks 

The power of the methods of theory of dynamical systems is due to their generality. The predictions we have 
made about the behavior of the field in the cavity are applicable to any motion of the mirror, not only to 
particular examples. Within our approach, we gave a complete classification of the possible resonances (phase 
locking) in the system, predicted that, generically, small detuning does not destroy the resonance, gave a 
simple explanation of the squeezing of the wave packets, interpreted the origin of the different contributions 
to the energy density in the cavity. We would also like to emphasize that our technique is non-perturbative. 

In the case of resonance, the standard numerical methods for solving partial differential equations would 
be very difficult to apply because of the concentration of the field in narrow packets. The proposed method. 
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Figure 5. Rescaled "shape" of a wave packet at times 2 (dashed hne), 3 (dot-dashed Une), 4, 6, 8, 10, 12 
(solid lines) (see the text). 

however, relies on iterating one-dimensional maps, so that resonances do not present any additional difSculty. 
The computer programs used to produce the pictures in this paper took minutes to run on a PC. 

Our methodology easily generalizes (see our paper [Hi]) to the case of quasiperiodic motion of the mirror, 
and the case of two moving mirrors (studied previously by Ji et al ^3 , Dodonov ^ Dalvit and Mazzitelli 
|10| , Li and Li 43 ) . Recently, similar ideas from dynamical systems have been applied to the study of waves 
in a fluid in two-dimensional basin by Manders et al j48| . 

Interestingly, the behavior of the field of the cavity (described by a partial differential equation) is easier 
to analyze than the behavior of a particle bouncing back and forth between two perfectly reflecting walls 
(assuming that the reflections are perfectly elastic). The latter system, suggested by Fermi |57j as a possible 
mechanism for acceleration of the particles in the cosmic rays, reveals a much richer dynamical behavior 
(see, e.g., the book of Lichtenberg and Lieberman :44; ). 

There are many questions that deserve a further study. An interesting question is whether theory of 
dynamical systems can be applied to the case of a constant-length cavity filled with dielectric with changing 
properties. Another problem is the absence of resonances of certain type (noticed in our paper |E] and 
discussed by W§grzyn |57| 'l - is it generic, and which resonances are forbidden? Can methods of dynamical 
systems be used to study the problem in higher dimensions? Can our methodology be applied to the recent 
suggestion by Jaffe and Scardicchio to apply methods of geometric optics to the study of Casimir effect? 
Can similar methods be applied to other fields (see, e.g., the study of a classical massive field in a pulsating 
resonator by Dittrich and Duclos [T?>])? 

The appearance of the Schwarzian derivative hints at possible deeper connections between the quantum 
problem of a moving mirror, partial differential equations, and dynamical systems. 
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